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We report on the calculation of the next-to-leading order QCD corrections to the production 
of top-antitop-quark pairs in association with a hard jet at the Tevatron and at the LHC. Results 
for integrated and differential cross sections are presented. We find a significant reduction of the 
scale dependence. In most cases the corrections are below 20% indicating that the perturbative 
expansion is well under control. Moreover, the forward-backward charge asymmetry of the 
top-quark, which is analyzed at the Tevatron, is studied at next-to-leading order. We find large 
corrections suggesting that the definition of the observable has to be refined. 
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1 Introduction 



The top-quark is by far the heaviest elementary fermion in the Standard Model (SM). With a 
mass of (172.6 + 1 .4) GeV [1] its mass is about 36 times larger than the mass of the next heaviest 
fermion, the bottom quark. The large mass has lead to various speculations whether the top- 
quark behaves as a normal quark or whether it plays a special role in particle physics. The 
electroweak SU(2)xU(l) gauge structure of the SM, which is successful in describing a large 
variety of measurements, would require the quark masses to be zero if the symmetry was un- 
broken. With the largest mass amongst the quarks it is thus natural to assume that the top-quark 
is most sensitive to the mechanism of electroweak symmetry breaking. In particular, the fact 
that the top-quark mass is close to the scale of electroweak symmetry breaking — or equivalently, 
that the Yukawa coupling to the Higgs is very close to one — has motivated different scenarios 
in which the top-quark drives the electroweak symmetry breaking. More details can be found 
in recent review articles [2, 3]. 

Ignoring the SM as the theory of particle physics one might still wonder whether the top- 
quark, which is almost as heavy as a gold atom, behaves as a point-like particle. A deviation 
from the point-like nature would appear as anomalous moments yielding differential distribu- 
tions different from the point-like case. Anomalous couplings to the gluon are most naturally 
probed via the production of an additional jet. An indirect measurement through the measure- 
ment of the total cross section is in general more difficult. This is in particular true when the 
interference term with the corresponding Born amplitude gives no contribution to the total cross 
section due to discrete symmetries. 

In the context of the SM we have the remarkable fact that the electroweak top-quark in- 
teractions are completely determined through the aforementioned SU(2)xU(l) gauge structure 
of the SM. The only free parameter appearing in top-quark physics is thus the top-quark mass 
or equivalently the Yukawa coupling to the Higgs boson. Once this parameter is measured all 
remaining properties are predicted. 

An important task for the ongoing Tevatron collider and the recently started LHC is the 
precise measurement of the top-quark properties. The ultimate goal is to measure the spin 
and the quantum numbers of the top-quark as precisely as possible. Any deviation from the 
SM would signal new physics. There is a variety of measurements which are currently done 
at the Tevatron and will be continued at the LHC. The total cross section, which is currently 
measured at the Tevatron with an accuracy of about 10% [4], is expected to be measured at the 
LHC with an accuracy of 5%. This measurement allows to make precise tests of the production 
mechanism. Other important measurements comprise the cross section for single top-quark 
production [5-11], the W-polarisation in top-quark decay or the spin correlations of top-quark 
pairs [12-20]. Of great interest is also the electric charge of the top-quark and its couplings to 
the Z-boson and the Higgs boson. They can be constrained via the measurements of the cross 
sections for tty, ttZ [21], and tt// [22-25] production. 

The production of a top-quark pair together with an additional jet is a further important 
reaction. This becomes already clear from the simple observation that a substantial number of 
events in the inclusive top-quark sample is accompanied by an additional jet. Depending on 
the energy of the additional jet the fraction of events with an additional jet can easily be of the 
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order of 10-30% or even more. For example at the LHC we find a cross section of 376 pb for 
the production of a top-antitop-quark pair with an additional jet with a transverse momentum 
above 50 GeV. This is almost half of the total top-quark pair cross section which is 806 pb [26] if 
evaluated in next-to-leading order For a more precise understanding of the topology of 

top-quark events it is thus important to have also an improved understanding of top-quark pair 
production together with a jet. As mentioned already above, this reaction provides a sensitive 
tool to search for anomalous top-quark-gluon couplings. The emission of an additional gluon 
also leads to a rather interesting property of the cross section: The differential cross section 
contains contributions from the interference of C-odd and C-even parts of the amplitude [27- 
30], where C denotes the charge conjugation (for a similar effect in QED see for example 
Refs. [31,32]). While for the total cross section these contributions cancel when integrating 
over the (symmetric) phase space they can lead to a forward-backward charge asymmetry of 
the top-quark which is measured at the Tevatron [33,34]. It should be stressed that no parity- 
violating interactions are involved. Note that the naively defined forward-backward charge 
asymmetry is zero at the LHC due to the symmetric initial state. A definition that leads to a 
non-trivial prediction here requires to select a preferred axis for each event [35], but it is not yet 
clear whether an asymmetry survives that is significant over all uncertainties. In inclusive top- 
quark pair production at the Tevatron the charge asymmetry appears first at one loop, because 
it results from interferences of C-odd with C-even parts of double-gluon exchange between 
initial and final states. The asymmetry for the inclusive sample has been studied in detail 
in Refs. [27-30]. The available predictions for tt production — although of one-loop order — 
describes this asymmetry only at leading-order (LO) accuracy in QCD. Recently the analysis 
has been extended to take large threshold logarithms at the next-to-leading-log (NLL) level into 
account [36]. It was found that at least this class of higher-order contributions do not change 
the theoretical prediction dramatically. The main reason is that this type of corrections affect 
the asymmetric cross section roughly in the same way as the symmetric one. In the ratio the 
corrections thus cancel to a large extent and lead to a stable theoretical prediction. In tt+jet 
production the asymmetry appears already in the tree amplitude. Thus, the NLO calculation 
described in this article provides a true NLO prediction for the asymmetry. The calculation 
presented in this work is an important tool in the experimental analysis of this observable at the 
Tevatron where the asymmetry is measured [33,34]. In a previous letter [37] we reported that 
the asymmetry receives large corrections. In this paper we study the situation in more detail for 
various values of the lower cut on the transverse momentum of the hard tagging jet. 

As aforementioned it is expected that the total cross section for top-quark pair production 
will be measured at the LHC with an accuracy of the order of 5%. Recently it has been shown 
in Refs. [26,38,39] that the accuracy of the currently available NLO predictions is only at 
the level of 12% (at NLO, but further reduced by the inclusion of the threshold logarithms) 
and largely dominated by the scale uncertainty. In Ref. [26] an estimate to the next-to-next- 
to-leading order (NNLO) cross section has been given. The approximation is based on the 

'Both numbers correspond to a top mass of 174 GeV. The number for top-antitop-quark pair production with 
an additional jet was obtained using CTEQ6M as pdf set, the one for the total top-quark pair cross section with 
CTEQ6.5 
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assumption that the NNLO corrections will be dominated by the threshold region as it is the 
case for the NLO corrections. In the threshold region the logarithmic behaviour together with 
the two-loop Coulomb singularity is derived from general arguments. In addition the complete 
scale dependence at two loops is included in the approximation. Using this approximation to the 
full NNLO result it is shown in Ref. [26] that the theoretical uncertainty may decrease to a few 
per cent. The remaining scale uncertainty is of the same order as the uncertainty induced by the 
parton distribution functions. Recently some progress towards a complete NNLO calculation 
has been made [4(M-8]. The one-loop corrections to tt + 1-jet constitute an important ingredient 
to the NNLO calculation of it production at hadron colliders. In this context we mention that — 
besides our NLO calculation presented here and in Ref. [37] — part of the one-loop amplitudes 
to gg — > ttg have also been evaluated in Ref. [49]. 

Apart from its significance as signal process it turns out that tt+ 1-jet production is also 
an important background to various new physics searches. A prominent example is Higgs 
production via vector-boson fusion. This reaction represents an important discovery channel for 
a SM Higgs boson with a mass of up to several 100 GeV [50,51]. The major background to this 
reaction is due to tt-l- 1-jet [52], again underlining the need for precise theoretical predictions 
for this process. 

It is well known that predictions at LO in the coupling constant of QCD are plagued by 
large uncertainties. In many cases the LO predictions in QCD give only a rough estimate. 
Only by including NLO corrections a quantitatively reliable prediction can be obtained. Given 
that the conceptual problems in such calculations are solved since quite some time, one might 
think that doing the required calculations should be a straightforward task. Unfortunately it 
turns out that this is not the case. The calculation of radiative corrections for 2 — > 3 and 2 — > 4 
reactions is still non-trivial — not speaking about reactions with an even higher multiplicity]^ 
The complexity of the corresponding matrix elements renders computer codes quite lengthy 
and CPU time consuming. The (more or less) automatically generated code may in addition 
lead to numerical instabilities. In particular, the reduction of one-loop tensor integrals to scalar 
one-loop integrals is in general difficult to do in a numerically stable way. In that context the 
calculation of the one-loop corrections to top-quark pair production with an additional jet is 
also interesting as a benchmark process for the development of new methods. 

In this paper we extend our previous work [37] on the NLO QCD corrections to tt+jet 
production at hadron colliders, where we discussed the scale dependence of the integrated cross 
sections at the Tevatron and the LHC and of the top-quark charge asymmetry at the Tevatron. 
We supplement this discussion upon including more numerical results showing the dependence 
on the lower cut set on the transverse momentum of the hard tagging jet and present first results 
on differential distributions. Moreover, we provide numerical results on the virtual one-loop and 
real-emission corrections for single phase-space points, in order to facilitate future comparisons 
to our calculation. 

This article is organized as follows: In Section [2] we briefly describe the calculation of the 
NLO corrections. Numerical results are presented in Section |3] In the appendices we provide 

2 More details and references on problems and suggested solutions can, for instance, be found in reports like 
Refs. [53,54]. 



3 



Figure 2.1: Representative sets of LO diagrams for gg fusion and qq annihilation in hadronic 
tt+jet production. 

numerical results on the virtual and real corrections for individual phase-space points; more- 
over, we collect the tables with the results for the differential cross sections there. 

2 Outline of the calculation 

2.1 Born approximation 

In Born approximation the partonic reactions are 

gg-»tfg, W-»ttg, ?g->% gq^iiq. (2.1) 

The last three reactions are related by crossing. Therefore, the required generic matrix elements 
are 

-> ttggg, -» ttqqg. (2.2) 

In the following we generically denote the external momenta and helicities with {/?,•} and {At} and 
identify the light partons with i = 1,2,3. The letter of a specific parton if used as an argument 
denotes the combination of spin, momentum, and (if relevant) colour of this parton, i.e. t = 
(p t ,A t ,i t ) or gi = (pi,Ai,ai). Representative sets of Born diagrams for the gg and qq channels are 
depicted in Figure [2T1 In total, there are 16 LO diagrams for — > ttggg and 5 for — > ttqqg. 
The colour decomposition for a tree amplitude corresponding to the process — > ttggg is 

tff(t,gi,g 2 ,g3,t) = ^(r^ir^r^^Afag^.g^.g^,!), (2.3) 

0-6S3 

where g s is the strong coupling constant, 53 the symmetric group and a = (0-1,0-2,0-2) e S3 
a permutation. The generators of the SU(N) gauge group in the fundamental representation 
are given by T a , and a, is the colour index of gluon g,-. The sum extends over all permuta- 
tions in 53. Physically this corresponds to all possible colour orderings of the gluons. The 
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function (t,gi,g2,g3,t) is thus the colour-ordered subamplitude often also called partial 
amplitude. Due to the colour-ordering, only diagrams with a particular ordering contribute to 
Ag 0) (t,gi,g2,g3,t). The partial amplitudes contain the kinematic information and are individu- 
ally gauge invariant. The corresponding colour decomposition for the process — > ttqqg reads 



ftf(\Xq,q,g) = gl 



1 



1 



(0), 



(2.4) 



We note that the amplitudes Al ? (i = 1,2,3,4) are linearly dependent. The relation 

U_A 5,l +/i 5,2 A 5,3 ^5,4 



can be used to express for example A^ in terms of A^ ( . with j = 1,2,3. This is particularly 
useful for the evaluation of the squared amplitude. Compact analytic results for the LO am- 
plitudes are given in Ref. [20] where the amplitudes have been used in the calculation of the 
NLO corrections for top-quark pair production. In addition we also performed several indepen- 
dent calculations, including one with Madgraph [55], and found complete numerical agreement 
among all those calculations. 



2.2 Virtual corrections 

The virtual corrections consist of the one-loop corrections to the LO reactions. One can classify 
the corrections into self-energy, vertex, box-type, and pentagon-type corrections where all the 
external legs are directly connected to the loop thus forming a pentagon. The latter are the most 
complicated ones due to their complexity and the involved tensor integrals. Typical examples of 
the pentagon graphs are shown in Figure [2T2l Specifically, there are 24 pentagons for — > ttggg 
and 8 for — > ttqqg. The total number of diagrams is 354 for the — > ttggg case and 94 for the 
— > ttqqg case. The challenging step in this context is the numerically fast and stable reduction 
of the tensor integrals to scalar one-loop integrals. 

Before describing the details we briefly outline the general setup. Owing to the involved 
kinematics the individual Feynman diagrams lead to large expressions which are cumbersome 
to evaluate. To be able to handle the large expressions and to ensure a fast numerical evaluation 
at the end we used a decomposition of the amplitude according to the spin and colour structure. 
Schematically the decomposition of the one-loop amplitude Jv^' reads 

= J]c c xSs({piU<ii})xf£\{pi-Pj}), (2.5) 

c,s 

where C c denote the colour structures, S s are the spin structures (elsewhere called "standard 
matrix elements"), and the functions are scalar functions that depend only on the scalar 
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Figure 2.2: Representative sets of pentagon diagrams for gg fusion and qq annihilation in 
hadronic tt+jet production at NLO QCD. 

products of the external momenta In detail, for — » ttggg there are 10 independent colour 
structures C c of the form 



' r 3 



if 



'hJV 



where the structure constant of the SU(3) gauge group f abc is defined in the usual way through 



[T a ,T°] =if abc T 



(2.6) 



We note that the counting refers to the group SU(3); for SU(N) there are 1 1 independent struc- 
tures. For N=3 the generic structure Tr[T ai T a2 T a3 ]6i t j t of the SU(N) case can be further re- 
duced, because d aia2ai 6 hh appearing through Tr[T° l T a2 T ai ] = l(d aia2 ° 3 + if aia2 ° 3 ) is express- 
ible in terms of the other 10 structures within SU(3) representations. Allowing only the use of 
relations which are compatible with conventional dimensional regularisation when simplifying 
the Lorentz structure of the amplitude, we find a few hundred spin structures, such as 

[v- t u t ](si -s 2 )(£3 -Pi), [f^3Wt](£i -S2), [i>t&&«t] Oi • P3), etc., 

where an obvious notation for the Dirac spinors Tr t , u t and gluon polarization vectors e,- is used. 
Restricting to four dimensions and using explicitly four-dimensional helicity techniques the 
number of independent structures may be further reduced. Each gluon polarization vector can 
be "gauged" to be orthogonal to an arbitrary light-like reference vector, a fact that reduces the 
algebraic expressions considerably. A useful choice is, e.g., given by the cyclic set of conditions 



P2-£\ =P3-S2=Pl-£3=0, 



(2.7) 



which supplements the trans versality relations pt ■ Si = 0. However, more than 100 spin struc- 
tures still remain (in D dimensions) in spite of this simplification. For — > tt^^g there are only 
4 colour structures (as in LO), 



C T^fl C r J~'Cl C 'Y'Cl 

°kh l i q ji' °'Wf i ids' °^h 1 i a 



hjq' 



tqjq 



(2.8) 
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but also more than 100 spin structures, such as 

[VtUt\[VqPtUqKe-p t ), [VtU t ] [Vq^Uq], [Vtf "t] [Vq7llUq](s • p t ), etc. . 

We stress that the algebraic reduction of each Feynman diagram to the standard form shown 
in Eq. (12.51) proceeds in D = 4 - 2e space-time dimensions, i.e. only Dirac equations, transver- 
sality and gauge conditions of polarization vectors, and momentum conservation are used. 
We have followed two different strategies for the evaluation of the decomposition shown in 
Eq. (12.51 ) . In one implementation all Feynman diagrams were first combined and then pro- 
jected onto the different structures. The idea behind this is to form gauge independent quantities 
where some cancellations may happen. In the second approach the decomposition is applied 
to each individual Feynman diagram. In the latter the numerical evaluation of the functions f cs 
is most efficiently done if the colour structures completely factorizes. For diagrams without 
4-gluon vertices this is trivially the case. Diagrams where the number of 4-gluon vertices 

is greater than are decomposed into 3" 4 terms, each with its own colour structure. Denoting 

CD 

such one-loop (sub)graphs generically T, their contributions / CiS to the functions f cs are written 

as 

fcV = c<PfP (2.9) 

with constants cy* , i.e. the total colour structure of Y is J^ c C c Cc- The colour-stripped functions 

CD 

fs , which contain the time- 
Writing the LO amplitude as 



CD 

fs , which contain the time-consuming loop functions are the same for all colour channels 



^ = J^CcA<® (2-10) 

c 

the contribution of Y to the one-loop-corrected spin- and colour-averaged squared amplitude is 
evaluated as follows, 

colour, spin 

= 2Re -^Ccccf^fPiiprpj}) ^A^({p ; },R}r 5,({p ; },R}) . (2.11) 

[Ai) ' J 

In detail, the colour correlation matrix 

C c& =YjC\C & (2.12) 

colour 

(0) * 

is calculated only once and for all for the whole process, and the interferences A 5 c S s of the LO 
amplitude and the different spin structures S s are only calculated once per phase- space point for 
every spin state. The approach where first all the diagrams are combined and then the projection 
is done yields similar formulae. The final formula then reads: 



£ 2Re{(tff)t* 5 



colour, spin 

= 2Re 



Z Z Z c ccfcs(iPi ■ Pj})ff s *dPi ■ Pj)) S s {{Pi), {Ai})S* s ,({ Pi }, {At}) \ , (2.13) 



j/t/j c,& s,s' 
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where / c y are the scalar functions appearing in the decomposition of the Born amplitude. To 
ensure the correctness of our results the two slightly different approaches were implemented in 
two complete independent computer codes. We note at this point that no significant difference 
between the two approaches concerning speed and numerical stability was observed. We also 
note that we used as far as possible different methods and also different tools to obtain the 
various ingredients discussed above. In the following we give some details of the techniques 
employed in the two implementations. 

Version 1 of the virtual corrections is essentially obtained following the method described 
in Ref. [23], where ttH production at hadron colliders was considered. Feynman diagrams and 
amplitudes have been generated with the FeynArts package [56, 57] and further processed with 
in-house Mathematica routines, which automatically create an output in Fortran. The infrared 
(IR), i.e. soft and collinear, singularities — which are treated in dimensional regularisation in 
both calculations — are analytically separated from the finite remainder in terms of triangle sub- 
diagrams, as described in Refs. [23,58]. This separation, in particular, allows for a transparent 
evaluation of so-called rational terms that originate from D-dependent terms multiplying IR di- 
vergences, which appear as single or double poles in (D-4). As generally shown in Ref. [59], 
after properly separating IR from ultraviolet (UV) divergences such rational terms originating 
from IR divergences completely cancel; this general result is confirmed in our explicit calcu- 
lation. The tensor integrals appearing in the pentagon diagrams are directly reduced to box 
integrals following Ref. [60]. (Similar methods have been proposed in Ref. [61].) This method 
does not introduce inverse Gram determinants in this step, thereby avoiding notorious numerical 
instabilities in regions where these determinants become small. Box and lower-point integrals 
are reduced a la Passarino-Veltman [62] to scalar integrals, which are either calculated analyt- 
ically or using the results of Refs. [63-65]. Sufficient numerical stability is already achieved 
in this way. Nevertheless the integral evaluation is currently further refined by employing the 
more sophisticated methods described in Ref. [66] in order to numerically stabilize the tensor 
integrals in exceptional phase- space regions. 

Version 2 of the evaluation of loop diagrams starts with the generation of diagrams and am- 
plitudes via QGRAF [67], which are then further manipulated with Form [68] and Maple and 
eventually automatically translated into C++ code. The reduction of the the 5-point tensor in- 
tegrals to scalar integrals is performed with an extension of the method described in Ref. [69]. 
In this procedure also inverse Gram determinants of four four-momenta are avoided. The 
lower-point tensor integrals are reduced using an independent implementation of the Passarino- 
Veltman procedure. The IR- finite scalar integrals are evaluated using the FF package [70,71]. 
Although the entire procedure is sufficiently stable, further numerical stabilization of the tensor 
reduction is planned following the expansion techniques suggested in Ref. [72] for exceptional 
phase-space regions. 

As stated above we used dimensional regularisation to regularise UV as well as soft and 
collinear divergences. We renormalised the coupling in a mixed scheme where the light flavours 
are treated according to the modified minimal subtraction MS , while the top-quark loop of the 
gluon self-energy is subtracted at zero momentum. The top-quark mass is renormalised in the 
on-shell scheme. More specifially we used the renormalisation constants as given for example 
in Ref. [23]. In these formulae the divergences of UV and IR origin are separated. They allow 
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us to check UV and IR finiteness separately. 



2.3 Real corrections 

The generic matrix elements for the real corrections are given by 

0->ttgggg, 0->tiqqgg, Q-> tlqqq'q', 0->tiqqqq (2.14) 

with q± q' . The various partonic processes are obtained from these matrix elements by all pos- 
sible crossings of light particles into the initial state. While the crossing symmetry is extremely 
helpful in constructing the required amplitudes it should be kept in mind that the large number 
of possible channels obtained from the different crossings lead to a significant increase in the 
computational complexity, given that every channel has to be integrated over the phase space. 
The amplitude for the process — » ttqqqq with identical quarks q can be obtained from the 
amplitude of the process — » tlqqq'q' with non-identical quarks q and q'\ 

&f\t,t, q,q,q,q) = ftf\t,l,q,q,q',q')-3lf\t,t,q,q',q',q). (2.15) 
The colour decomposition of a tree amplitude corresponding to the process — > ttgggg is 

4 0) (t,gi,g 2 ,g3,g4,t) = ^2 ( r ^ 7 ^ r ^ r<fe ' 4 ) Ut A2 )) (t,g . 1 ,g^,g^,g ff . 4 ,^ (2.16) 

cre5 , 4 

where the sum is over all permutations cr = (0-1,0-2,0-2,0-4) of the symmetric group 54. The 
colour decomposition for the process — > tlqqgg reads 

4 r 

^ 0) (t,t,^, gl ,g 2 ) = I J] [*U|(^^)^©t.l*.«.^.g^) 

(a- 1 a-2)eS2 

+r 2 T lk A Z ia q > q > ' + ^ T ° LT2 Xk 6 vAZ^ q > q > &T, . foi) 



_ I (7*1 7^2 ) . . 5. Ag(t, t, <?, g^ , g a2 ) 



(2.17) 



where again the sum is over all permutations of the gluon legs. Finally, the colour decomposi- 
tion of the process — > tlqqq'q' is 

Mf\t,t,q,q,q',q') = -j 8 i j-6 iqjq ,6 iq>j - t A ( ^ l (t,l,q,q,q' ,q') + 6 kjq ,6 iq ,j q 6 iqj - t A^ l (t,l,q , ,q , ,q,q) 

-^u Aa<V ? ' a 6 2^' l > 9, q,q',q')- ^ kj . ql s iq/ ^^(t, t, q ',q', q ,q) 
-^ 6 ^k/jAji A Z^ q ^^^ q,) + ^ 6 hrA,jA'^ A Z {t ~ Uq ^ q ^ q ' q) ■ (2 - 18) 
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To extract the IR singularities and for their combination with the virtual corrections we 
employ the dipole subtraction formalism [73-75] . Specifically, the formulation [75] for massive 
quarks is used. At NLO schematically one has the following contributions: 



(O) 



NLO 



= J O n+ \dcr R + J O n dcr v + J O n dcr c . 



(2.19) 



77+1 



Here dcr R denotes the real emission contribution, whose matrix elements are given by the square 
of the Born amplitudes with 6 partons |^g 0) | 2 , dcr v is the virtual contribution, whose matrix 

elements are given by the interference term of the one-loop amplitudes jTll 1 ^ with 5 partons 

with the corresponding Born amplitude and dcr c denotes a collinear subtraction term, 
which originates from the factorisation of the initial-state collinear singularities. The function 
O n defined on the n-particle phase space stands for any prescription (^-functions for phase- 
space cuts, ^-functions for distributions) defining an IR-safe observable. Taken separately, the 
individual contributions are IR divergent, and only their sum is finite. In order to render the 
individual contributions finite, so that the phase- space integrations can be performed by Monte 
Carlo methods, one adds and subtracts a suitably chosen "counterterm" dcr A : 



(O) 



NLO 



O n+ id(T R -O n do-' 



n+1 



O n d(T V + O n d(T C + O n J 

1 



da 1 



(2.20) 



The matrix element corresponding to the approximation term dcr A is given as a sum over 
dipoles: 



pairs i,j ktij 



(2.21) 



Each dipole contribution has the following form 
1 



T* T 



^ ) *(p\,...,p(ifr...,pk,.-)—^Vij > k&f 



T 2 

'j 



Here T, denotes the colour charge operator for parton i and Vuk is a matrix in the spin space 
of the emitter parton (z'j). The momenta p^f) and pk are obtained from the momenta puPj 
and pk. In general, the operators T, lead to colour correlations, while the Vy^'s may lead to 
spin correlations. The approximation dcr A has to fulfill the requirement that dcr A is a proper 
approximation of dcr R with the same point- wise singular behaviour (in D = 4 - 2e dimensions) 
as dcr R itself. Thus, dcr A acts as a local counterterm for dcr R , and one can safely perform the 
limit s —> 0. This defines the finite contribution 



(2.22) 



n+1 
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The subtraction term can be integrated over the unresolved one-parton phase space. Due to this 
integration, all spin correlations average out, but colour correlations still remain. In a compact 
notation, the result of this integration is often written as 

dcr c + J dcr A = I<g>dcr s + K® da B + V® da B . (2.23) 
1 

The notation ® indicates that colour correlations still remain and that an integration is involved. 
The term I ® da B lives on the phase space of the Born configuration and has the appropriate 
singularity structure to cancel the IR divergences coming from the one-loop amplitude. There- 
fore, dcr v + I®dcr B is IR finite. The terms (K + P) <8> dcr B involve in addition an integration over 
the momentum fraction x that rules the collinear splitting of the incoming parton. From the 
integration of the subtraction terms we obtain the finite contribution 



J O n (do y + \®do- B + K®d(T B + V®d(T B ) E=Q . (2.24) 



The final structure of an NLO calculation in the subtraction formalism is then 

<0> NLO = (0)™f 1] + (0)^°. (2.25) 

Since both contributions on the right-hand side of Eq. (12.251) are now finite, they can be eval- 
uated with numerical methods. The explicit forms of the dipole terms Duk, together with 
the integrated counterparts, can be found in Ref. [73] for massless QCD and in Refs. [74,75] 
including massive quarks. 

Analogously to our evaluation of the virtual corrections, we have also performed two inde- 
pendent calculations of the real corrections. 

One calculation of the real corrections results from a fully automated calculation based 
on helicity amplitudes, as described in Ref. [76]. Individual helicity amplitudes are computed 
with the help of Berends-Giele recurrence relations [77]. The evaluation of colour factors and 
the generation of subtraction terms is automated. For the channel gg — » ttgg a dedicated soft- 
insertion routine [78] is used for the generation of the phase space. 

The second calculation uses for the LO 2 — > 3 processes and the gg — » ttgg process optimized 
code obtained from a Feynman diagrammatic approach. As in the calculation described before, 
standard techniques like colour decomposition and the use of helicity amplitudes are employed. 
For the 2 — > 4 processes including light quarks, Madgraph [55] has been used. The subtraction 
terms according to Ref. [75] are obtained in a semi-automatized manner based on a library 
written in C++. 

The two independent computer codes were compared point- wise at a few phase-space points. 
In addition the entire numerical integration of the real corrections was done independently us- 
ing the two codes. We found complete agreement of the numerical results when the numerical 
uncertainty from the phase-space integration is taken into account. 
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3 Numerical results 



3.1 Setup 

In the following we consistently use the CTEQ6 [79] set of parton distribution functions (PDFs). 
In detail, we take CTEQ6L1 PDFs with a one-loop running a s in LO and CTEQ6M PDFs with 
a two-loop running a s in NLO. The number of active flavours is Afc = 5, and the respective QCD 
parameters are A^ = 165MeV and A^ s = 226 MeV. As mentioned earlier the top-quark loop 
in the gluon self-energy is subtracted at zero momentum. In this scheme the running of a s is 
generated solely by the contributions of the light quark and gluon loops. The top-quark mass 
is renormalized in the on-shell scheme, as numerical value we take m t = 174 GeV. If not stated 
otherwise, we identify the renormalization and factorization scales, yu ren and //f aC b with m t . 

For the definition of the tagged hard jet we apply the jet algorithm of Ref. [80] with R = 1 
and require a transverse momentum of 7?T,jet > PT,jet,cut with PT,jet,cut = 20 GeV and PT,jet,cut = 
50 GeV for the hardest jet at the Tevatron and the LHC, respectively. The outgoing (anti)top- 
quarks are neither affected by the jet algorithm nor by the phase-space cut. We assume them as 
always tagged. Note that the LO prediction and the virtual corrections are not influenced by the 
recombination procedure of the jet algorithm, but the real corrections are. 

Up to the transverse-momentum cut /?T,jet,cut = 50 GeV for the LHC, the setup used in this 
article coincides with the one used in Ref. [37]. There, PT,jet,cut = 20 GeV was used both for the 
Tevatron and the LHC. 

3.2 Results for the Tevatron 

As discussed in Ref. [37], the integrated LO cross section for tt-l- 1-jet production at the Teva- 
tron is dominated by the qq channel with about 85%, followed by the gg channel with about 7%. 
This is rather similar to inclusive top-quark pair cross section where again at LO about 90% is 
obtained from qq channel and about 10% from the gg channel. In difference to the inclusive 
case the qg (qg) are not suppressed in the coupling. This accounts for the slightly larger contri- 
bution from these channels. One should keep in mind that the precise contribution of individual 
channels depends on the factorisation scale as well as on the chosen parton distributions. As a 
consequence the aforementioned numbers give just a qualitative picture. 

In Table 13. ll we provide the LO and NLO predictions for the integrated cross sections for 
different values of the cut on the transverse momentum of the hard jet (left part). The values 
presented are for the central scale n = yUf act = // ren = wi t . In parentheses we quote the uncertainty 
due to the numerical integration. The scale dependence is indicated by the upper and lower 
indices. The upper (lower) index represents the change when the scale is shifted towards n = 
m t /2 (// = 2m t ). Rescaling the common scale /u = yUf act = /i re n from the default value m t up (down) 
by a factor 2 changes the cross section in LO and NLO by about 60% (35%) and 9% (18%), 
respecively, i.e. the scale uncertainty is reduced considerably through the inclusion of the NLO 
corrections. The above findings are rather insensitive to the chosen cut value. We find only 
variations at the per-cent level. In particular, there is no big difference for the lowest cut value 
compared to the other values, suggesting that the perturbative expansion is under control and 
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PT,jet,cut L^cvJ 


°"ttje 
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ttpb] 

NT n 


t n 
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nt n 

IN LU 


20 
30 
40 
50 


1.583(2)+°;* 
0.984(l)+°;§° 
0.6632(8)+°,;^ 
0.4670(6)^29 


1.791(l)+gif 
1.1194(8)+°;iJ 
0.7504(5)!«-^ 2 
0.5244(4)^-049 


-7.69(4)tJ0 5 
-8.29(5)+°g 5 
-8.72(5)!g;|3 
-8.96(5)!j; 4 


-1.77(5)*$ 
-2.27(4)+°;3| 
-2.73(4)^-35 
-3.05(4)^-49 



Table 3.1: Cross section cr t ij el and forward-backward charge asymmetry Ap B at the Tevatron for 
different values of />T,jet,cut for /j. = //f act = /z ren = m t . The upper and lower indices are the shifts 
towards jx = m t /2 and jx = 2m t . 



not spoiled by the appearance of large logarithms. Compared with the total cross section we 
find that for the small pj cut of 20 GeV the tt + 1-jet events represent almost 30 % of the total 
cross section. This fraction shows an evident dependence on the value chosen for the pj cut. 
The fraction is reduced to about 8% when 50 GeV is chosen for the cut. The NLO corrections 
change the ratio for a given value of pj cut only at the level of a few per cent. 

In the right part of Table l3.!l we show results for the forward-backward charge asymmetry. 
In LO the top-quark charge asymmetry is defined by 

^FB.LO = "Z+ ' (3.1) 

with the definition 

c^lo = ^Lo()>t>0) ± cr LO (y t <0), (3.2) 

where y t denotes the rapidity of the top-quark. Cross-section contributions cr(y t < 0) correspond 
to top-quarks in the forward or backward hemispheres, respectively, where incoming protons 
fly into the forward direction by definition. Denoting the corresponding NLO contributions to 
the cross sections by 5cr* L0 , we define the asymmetry at NLO by 

FB,NLO - a+ 1+ - + ' ^-3) 

LO \ LO LO ' 

i.e. via a consistent expansion in a s . Note, however, that the LO cross sections in Eq. (13.31) are 
evaluated in the NLO setup (PDFs, a s ). In Ref. [37] it was already pointed out that the LO asym- 
metry for a px-cut of 20 GeV, which is about -7.7% with a small scale uncertainty, is reduced 
to about -1.8% with the rather large scale uncertainty that is — assessed conservatively — almost 
as large as its absolute size. The reason for this growing scale uncertainty when going from LO 
to NLO simply results from the fact that the LO prediction for Ap B is independent of the renor- 
malization scale, since the strong coupling drops out in the ratio. Thus, the scale dependence 
does not reflect the total theoretical uncertainty in LO at all for this quantity. Table 13. II shows 
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that this feature qualitatively holds true also for larger values of pT,jet,cut use d for the jet defini- 
tion, but the LO and NLO asymmetries are shifted towards larger absolute values for a larger 
cut. 

Figure 13.11 shows the distributions in the transverse momenta of the hard jet, ;?T,jet> of the 
total ft system, /?x,tt> and of the top-quark, pj t . In LO PTjet an d PT,tt coincide because of mo- 
mentum conservation in the transverse plane, but the radiation of two jets in the real corrections 
renders them different. The numerical results in Figure [3TTT however, reveal that the differences 
are very small. The shown pj distributions drop with growing pj, where the spectrum for the 
top-quark is much harder than the ones for the jet and the tt system. In the PTjet and pj t i spectra 
93% of the events are concentrated below a pi of about 100 GeV at NLO, while 92% of the 
events have a px,t with less than 200 GeV. Employing a fixed scale p = m t , the NLO corrections 
do not simply rescale the LO shape, but induce distortions at the level of some 10%, which 
redistribute events from larger to smaller transverse momenta. We believe that two effects con- 
tribute to these distorsions. First of all, the use of a fixed renormalization scale p ren = m t is 
not an appropriate choice for high pj events. Due to the large value of a s the LO calculation 
overestimates the cross section for high pj and the NLO calculation has to compensate this 
scale choice. We expect that the distortions due to this effect are reduced if an appropriate 

/?T-dependent scale choice is used, such as p = ^jpj + mf. As a second effect in particular the 
pj distribution of the top-quarks can become softer due to the emission of additional particles, 
which is accounted for the first time by an NLO calculation. 

As can also be seen from the lower panel of each plot, we find again an important reduction 
of the scale dependence when the NLO corrections are taken into account. At least for the 
PTjer an d /?T,trdistribution the corrections are of moderate size. For the 7?T,t-distribution we 
find large corrections for large values of px,t- The corrections are almost 50% for a pj arround 
400 GeV. As mentioned this could probably be cured by employing a px-dependent scale. 

Figure I3T21 depicts the distributions in the pseudo-rapidity and rapidity of the top-quark, T] t 
and y t , and in the rapidity, yj et , of the hard jet. For massless momenta the pseudo-rapidity, which 
is defined through 

'0 



77 = -ln^-jj (3.4) 

(# is the scattering angle with respect to the beam axis), is equivalent to the rapidity defined 
through 

y-h4^\ (3.5) 



2 \E-p z 

(E denotes the energy, p z the three-momentum component along the beam axis). For the mas- 
sive top-quark we observe a rather important difference between the rj t and y t distribution. Re- 
cently there has been significant interest in the rapidity distribution of the jet, as MC@NLO [81] 
and Alpgen with MLM matching [82] disagree on this distribution. Our result includes for the 
first time the full 0(a A s ) matrix elements. Note that the set-up of ref. [82] differs from the one 
used here, therefore the distributions should not be compared directly. 

At NLO, 90% of the events are concentrated within |?7t| < 2.0. Demanding \y t \ < 1.2 selects 
96% of the events. For Vj et we find that 94% of the events are contained in [yj et | < 2.4. 
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Figure 3.1: Transverse-momentum distributions of the hard jet (pTjetX of the total tt sys- 
tem (px,tt)> and of the top-quark (pT,t) at the Tevatron. The lower panels show the ratios 
K = NLO/LO as well as the LO and NLO scale uncertainties corresponding to a rescaling 
of yu = //fact = jUren = m t by a factor 2. 
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Figure 3.2: Distributions in the pseudo-rapidity (rj t ) and rapidity (y t ) of the top-quark, and in the 
rapidity (jj et ) of the hard jet at the Tevatron. The lower panels show the ratios K = NLO/LO as 
well as the LO and NLO scale uncertainties corresponding to a rescaling of fi = fif act = /i ren = m t 
by a factor 2. 
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Table 3.2: Cross section cr t (j et at the LHC for different values of /?T,jet,cut for p = /if ac t = /^ren = mt- 
The upper and lower indices are the shifts towards p = m t /2 and p = 2m t . 



The reduction of the forward-backward asymmetry Ap B discussed above induced by the 
NLO corrections is clearly visible in the r/ t and y t distributions. The corrections are larger in 
the forward direction. The asymmetry in the LO distributions is thus reduced by the NLO 
corrections. It is hardly conceivable that this higher-order effect can be absorbed into LO 
predictions by phase- space-dependent scale choices. It should be realized that the forward- 
backward- symmetric rapidity distribution of the hard jet gets distorted by the corrections as 
well. The corrections increase for large values of |yj e t|. 

At least in the regions of the distributions in which the rate is not too much suppressed, 
the NLO corrections reduce the scale uncertainty of the LO distributions in a similar way as 
observed for the integrated cross section. 

3.3 Results for the LHC 

Table [3721 shows the integrated cross section for various values of the cut PT,jet,cut on the trans- 
verse momentum of the hard tagging jet at the LHC. In contrast to the Tevatron, the gg channel 
comprises about 70% of the LO pp cross section, followed by qg with about 22% [37]. We 
note that the importance of the qg channel is very different from the inclusive top-quark pro- 
duction. For inclusive top-quark pair production this channel is suppressed — despite the large 
parton luminosity in this channel — because it appears only at NLO. For tt+ 1-jet production 
the qg channel appears already in LO and thus gives a significant contribution due to the large 
parton luminosity. Comparing the LO and NLO predictions we find again that the large scale 
dependence of about 100% in the LO cross section is considerably reduced after including the 
NLO corrections. The ratio of the NLO tt+ 1-jet cross section to the total NLO tt cross section 
is about 47%, 22%, and 7% for a p T cut of 50 GeV, 100 GeV, and 200 GeV, respectively. 

In Figure [331 we show the distributions in the transverse momenta of the hard jet, pxjeb of 
the total tt system, p^ n, and of the top-quark, /?x,t- The distributions become harder in pj when 
going from the Tevatron to the LHC, as expected from the higher scattering energy. At NLO, 
92% of the events have transverse momenta PTjet < 250 GeV, and 94% have px,t < 300 GeV 
in the respective distributions. In contrast to the Tevatron, the PTjet and pj ;t j distributions, 
which are identical in LO become different in NLO. For the pTjerdistribution the lowest bin 
(0 < PTjet < 50 GeV) is always empty due to the cut applied. For the px,tt this holds also true 
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Figure 3.3: Transverse-momentum distributions of the hard jet (j?T,jetX of the total tt system 
(PT,tt)> and of the top-quark (px,t) at the LHC. The lower panels show the ratios K = NLO/LO as 
well as the LO and NLO scale uncertainties corresponding to a rescaling of yU = //f act = jx Kn = m t 
by a factor 2. 
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in LO because the transverse momenta between the tt system and the additional hard jet are 
balanced. In NLO the lowest bin in the pi a distribution is populated due to an additional jet. 
For large pj the difference between j^Tjet and /? Tt j distributions is at the level of about 10%. 
This is again due to the presence of the additonal jet. As observed already for the Tevatron, the 
shapes of the pi distributions receive distortions by the corrections. 

Figure I3T41 illustrates the distributions in the pseudo-rapidity and rapidity of the top-quark, 
77 t and y t , and in the rapidity, Vj et , of the hard jet. In the respective distributions 92% of the NLO 
events concentrate within the regions |^ t | < 3.0, 96% have |y t | < 2.4, and 96% have |yj et | < 3.6 
in the respective distributions, i.e. these distributions get broadened roughly by one unit in 
the transistion from Tevatron to LHC. As for the Tevatron, we find distortions of the shapes 
induced by the corrections that are hard to mimic by phase- space-dependent scale choices. All 
the shown y and rj distributions at the LHC are forward-backward-symmetric, but actually the 
distributions of top and antitop-quarks are intrinsically different. Numerically we do not observe 
a significant difference, so that we show only the distributions for the top-quark. 

Again, in all shown distributions a reduction of the scale uncertainty by the NLO corrections 
is visible that is comparable to the one in the integrated cross section. 

4 Conclusions 

The production of tt+jet final states represents important processes both at the Tevatron and the 
LHC. The signal is interesting in its own right, because large fractions of the it samples show 
additional jet activity and deviations from the SM could signal new physics such as top-quark 
compositeness. Moreover, tt+jet production delivers a large background to many searches at 
the LHC, such as for the Higgs boson via weak- vector-bo son fusion. 

We have presented NLO QCD predictions for tt+jet production at the Tevatron and the LHC. 
The NLO corrections reduce the scale uncertainty of the total cross section and of the differen- 
tial distributions compared to a LO calculation, which can only provide qualitative predictions. 
Further theoretical improvements could only be achieved by dedicated QCD resummations, 
since a full treatment at NNLO is certainly out of reach. Already the presented NLO calcula- 
tion is quite complicated. For this reason we have also documented a set of numerical results 
for the one-loop correction and for the real-emission parts at single phase-space points, in order 
to facilitate comparisons to our calculation by other groups. 

The charge asymmetry of the top-quark, which is measured at the Tevatron, is significantly 
decreased at NLO and is almost washed out by the residual scale dependence. We have studied 
the dependence of the NLO asymmetry on the cut on the transverse momentum of the hard 
tagging jet. Further refinements in the description of the charge asymmetry are required to 
stabilize the predictions with respect to higher-order corrections. Moreover, the top-quark decay 
should be taken into account properly. 

Finally, the presented NLO QCD calculation for the tt+jet process represents a building 
block for a full NNLO QCD prediction for tt production, a demanding and important calculation 
that is currently in progress by various groups. 
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Figure 3.4: Distributions in the pseudo-rapidity (rj t ) and rapidity (y t ) of the top-quark, and in 
the rapidity (yj e t) of the hard jet at the LHC. The lower panels show the ratios K = NLO/LO as 
well as the LO and NLO scale uncertainties corresponding to a rescaling of /u = /Uf act = // ren = m t 
by a factor 2. 
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Appendix 



A Benchmark numbers for the virtual corrections 

In order to facilitate a comparison to our calculation, we provide explicit numbers on the 
squared LO amplitude and the corresponding virtual corrections for a single non-exceptional 
phase-space point. The set of momenta for ab — > ttc with the explicit partonic reactions gg — > ttg, 
qq — > ttg, qg —> ttq, and gq — > ttq is chosen as 

p a = (500,0,0,500), 

Pb = (500,0,0,-500), 

p t = (458.5331753852783,207.0255169909440,0,370.2932732896167), 

Pl = (206.6000026080000,-10.65693677252589,42.52372780926147,-102.3998210421085), 

p c = (334.8668220067217,-196.3685802184181,-42.52372780926147,-267.8934522475083), (A.l) 

with the obvious notation p = (p° , p 1 , p 2 , p 3 ) and all the components given in GeV. The top- 
quark mass is set to m t = 174 GeV. We give numbers on the spin- and colour- averaged squared 

(0) ? (0) (1)* 

LO amplitude ffl 5 | as well as for the contribution 2Re(^ Jv 5 ). For the Born amplitude 
we factor out the coupling, we define 



The factor 1/Af c is due to the average over the incoming colour. For the channel gg, qq, qg, gq 
we have N c = 64, 9, 24, 24. Note that the coefficient ao only depends on the chosen phase-space 
point — which implicitly also contains the information about the top-quark mass. The results 
for ao are shown in Table lA.il Version 1 and Version 2 correspond to our implementations. For 
completeness we compare also with Madgraph. The one-loop contribution is renormalised 
and thus UV finite. However, the virtual corrections still contain collinear and soft singularities. 
Similar to what has been done in Ref . [54] in a comparison of the virtual corrections to WW+jet 
production, we use the decomposition 



The results are shown in Table I A. 21 In addition we also give the corresponding results for the 
I-operator of the dipole subtraction function as defined in Ref. [75], with the auxiliary param- 
eter k = 2/3. We us the same decomposition as for the one-loop corrections. The individual 
coefficients are shown in Table lA.3l Note that the coefficients c ; with i = -2,-1,0 contain one 
factor of ff s . We use 




(A.2) 




a s {p) = 



An 
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Table A. 1 : Colour and spin averaged LO matrix elements squared. 



with Nf = 5 and A = 226 MeV, leading to 

a s (m t ) = 0.1075205492734706. (A.5) 

For the LO amplitudes we find an agreement of at least 14 digits — pretty close to what one 
can get using 64bit double precision with 53bits for the mantissajj For the one-loop corrections 
we find at least an agreement of 10 digits for the finite terms. For the qq, qg, and gq channels, 
which are numerically less involved, we find an agreement of up to 14 digits. The finite terms 
from the I-operator agree even better. The coefficients of all divergences typically agree to 10 
digits, or better in the cases with external quarks. We note, however, that we do not cancel 
the IR divergences numerically. We also observe a cancellation between the finite part of the 
I-operator and the corresponding virtual contributions. For the gg — > ttg channel almost three 
digits are cancelled in the combination. 

3 Note that some of the code was run on a x86-64 architecture, where usually floating point arithmetic is typ- 
ically done in the SMD unit of the processor and thus restricted to 64 bit also for intermediate results. (On x86 
architectures arithmetic is usually done in the FPU which usually works with extended precision for intermediate 
results (80 bit). The precision is reduced when results are stored back to memory.). 
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Table A. 2: Coefficients for colour and spin averaged virtual corrections. 
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0.1540118421074569 


-0.0731096894943435 


-0.5280576886301999 


Version 2 


0.1540118421074573 


-0.0731096894943437 


-0.5280576886302015 


qq ttg 


Version 1 


0.0969704191046952 


0.0126983208241661 


-0.3992776407671517 


Version 2 


0.0969704191046950 


0.0126983208241662 


-0.3992776407671513 


qg ttq 


Version 1 


0.0969704191046950 


0.0056430956994063 


-0.4069645466913195 


Version 2 


0.0969704191046951 


0.0056430956994064 


-0.4069645466913194 


gq ttq 


Version 1 


0.0969704191046950 


-0.0833362739127882 


-0.3392937280293060 


Version 2 


0.0969704191046950 


-0.0833362739127883 


-0.3392937280293059 



Table A. 3: Coefficients for colour and spin averaged results for the I-operator. 
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B Benchmark numbers for the subtraction terms 



In this section we give also results for the subtraction term in the dipole formalism for one 
phase-space point. We find this useful to facilitate the comparison of upcoming calculations. 
Recently some effort has been invested to automatize this part of the calculation [83-86]. The 
results presented here may provide an interesting benchmark point for these attempts. The set 
of momenta for ab — » tied is chosen as 

p a = (2100,-0,-0,2100), 

Pb = (2800,-0,-0,-2800), 

p t = (1581. 1 18367308447, 1254.4623 16247655 -766.9360998604944,-554.7905976902205), 

p l = (1460.449317799282,-975.9731477430979,-466.5314749495881,965.6402060944737), 

p c = (545.4084744819,218.7220720302516,472.0439121434804,-163.7241712507502), 

p d = (1313.023840410371,-497.2112405348086,761.423662666602,-947.1254371535031), (B.l) 

with all components given in GeV. The top-quark mass is set to m t = 174 GeV. We give numbers 
on the spin- and colour-averaged squared real emission amplitude IJ?!^! 2 as well as for the sum 
of the subtraction terms. We define the numbers bo by 



The factor 1 /4 accounts for the average over the spins of the initial partons, the factor 1 /N c is 
due to the average over the colour of the incoming partons. N c contains a factor 8 for every 
incoming gluon and a factor 3 for every incoming quark or antiquark. S is the symmetry factor 
accounting for identical particles in the final state. The numbers for the dipole subtraction terms 
are defined analogously: 



pairs ij kj^ij 

Note that the factor for the average over the colour of the incoming partons is included in the 
definition of Dij,k- The numbers bo and do contain the strong coupling constant. As numeri- 
cal value for a s we use again Eq. (IA.4I) with Nf = 5 and A = 226 MeV, i.e. the value given in 
Eq. (IA.5I) . The results for bo and do are shown in Tables IBTI and IBT21 The two implementa- 
tions agree at least to 14 digits for the matrix elements squared and at least to 12 digits for the 
sum of the subtraction terms. 



In this appendix we give the tables for the differential distributions. For each distribution, we 
list the NLO predictions for the scale choice p. = mJ2, p = m t and p = 2m t . In all tables we have 
set p = p ien = /ifact- The errors result from the Monte Carlo integration. The bin is specified by 




spin, colour 



(B.2) 




(B.3) 



C Tables for histograms 
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MGeV- 4 ] 


JotGeV -4 ] 


g(Pa)g(Pb) -> t(Pt)t(Pi)g(Pc)g(Pd) 


Version 1 
Version 2 


3.12815868347843 • 10" 9 
3.12815868347842- 10" 9 


4.1037601540955 • 10" 9 
4.1037601540962- 10" 9 


q(Pa)q(Pb) -» t(Pt)t(Pt)g(Pc)g(Pd) 


Version 1 
Version 2 


4.48308845446477 • 10" 10 
4.48308845446475 • 10" 10 


4.90476067759639 • 10" 10 
4.90476067759631 • 10" 10 


q(Pa)g(Pb) -» t(Pt)t(Pt)g(Pc)q(Pd) 


Version 1 
Version 2 


1.10256509258713 • 10" 10 
1.10256509258713 • 10" 10 


1.919073353538 • 10" 10 

1.919073353539 • 10" 10 


q(Pa)g(Pb) -» i(PtKPt)q(Pc)g(Pd) 


Version 1 
Version 2 


1.384600673183816- 10" 10 
1.384600673183812- 10" 10 


3.3382231835799 • 10" 10 
3.3382231835798 • 10" 10 


g(Pa)g(Pb) -» t(PtNPt)q(Pc)q(Pd) 


Version 1 
Version 2 


2.42841040229558 • 10" 10 
2.42841040229557 • 10" 10 


4.271065781530- 10" 10 
4.271065781532- 10" 10 



Table B.l: Colour and spin averaged real emission matrix element squared and dipole subtrac- 
tion terms related to the processes — » ttgggg and —> ttqqgg. 



its central value. The bin width — which we chose constant for the entire histogram — is obtained 
from the distance of two neighboring bin positions. Note that we use the same definition for 
the cross section as described in Section [3] In particular, we demand a minimum pj for the 
additional jet. For the Tevatron 20GeV is used while for the LHC 50GeV is used. 
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^o[GeV- 4 ] 


J [GeV- 4 ] 


q(Pa)q'(pb) -> t(pt)t(pdq'(Pc)q(pd) 


Version 1 


4.44137855516180- 10" 12 


1.6381811832266 


10" 


ii 


Version 2 


4.44137855516180- 10" 12 


1.6381811832275 


10" 


■11 


q(Pa)q'(pb) -* t(pt)t(pi)q(pc)q'(pd) 


Version 1 


1.733763330485899- 10" 11 


1.06832579841007 


•10 


-10 


Version 2 


1.733763330485899- 10- 11 


1.06832579841000 


•10 


-10 


q(Pa)q'(pb) -> t(pt)t(pi)q'(Pc)q(pd) 


Version 1 


4.796260245409952 • 10~ 12 


1.8776008214791 


10" 


■11 


Version 2 


4.796260245409957 • 10~ 12 


1.8776008214799 


10" 


■li 


q(Pa)q(Pb) -* t(ptNpt)q'(Pc)q'(Pd) 


Version 1 


6.13924303047741 -10" 11 


6.990891152615- 


io- 


n 


Version 2 


6.13924303047739- 10- 11 


6.990891152614- 


10" 


n 


q(Pa)q(Pb) -> t(pt)t(Pt)q(Pc)q(Pd) 


Version 1 


1.371477814148721 • 10 -11 


4.1848434402744 


10" 


■li 


Version 2 


1.371477814148719- 10 -11 


4.1848434402750 


10" 


■li 


q(Pa)q(Pb) -* t(ptNpt)q(Pc)q(Pd) 


Version 1 


1.411042000289490- 10" 11 


6.3674516988854 


10" 


-li 


Version 2 


1.411042000289488- 10" 11 


6.3674516988857 


10" 


■li 


q(p a )q(Pb) -> t(ptNpt)q(Pc)q(Pd) 


Version 1 


2.054843839960259- 10" 11 


3.6253236096328 


10" 


■li 


Version 2 


2.054843839960252- 10" 11 


3.6253236096334 


lO- 


■li 



Table B.2: Colour and spin averaged real emission matrix element squared and dipole subtrac- 
tion terms related to the processes — > ttqqq'q' and — > ttqqqq. 
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Ptm [GeV] 


// = m t /2 


dcr r fb n 
^Tjet L GeV J 

H = m t 


yu = 2m t 


12.5 


16.83+0.07 


15.70 + 0.04 


13.4 + 0.2 


37.5 


38.47 + 0.04 


35.08 + 0.03 


28.83 + 0.02 


62.5 


12.67 + 0.02 


11.45 + 0.01 


9.340 + 0.009 


87.5 


5.29 + 0.01 


4.805 + 0.008 


3.92 + 0.005 


112.5 


2.465+0.007 


2.277 + 0.006 


1.864 + 0.003 


137.5 


1.216 ±0.005 


1.146 ±0.003 


0.943 + 0.002 


162.5 


0.629 + 0.004 


0.604 + 0.002 


0.496 + 0.001 


187.5 


0.326 + 0.003 


0.324 + 0.002 


0.2659 + 0.0007 


212.5 


0.173 + 0.002 


0.174 + 0.001 


0.1460 + 0.0005 


237.5 


0.093 + 0.001 


0.0945 + 0.0008 


0.0796 + 0.0004 


262.5 


0.047 + 0.001 


0.0522 + 0.0006 


0.0445 + 0.0003 


287.5 


0.0252 + 0.0007 


0.0285 + 0.0003 


0.0238 + 0.0002 


312.5 


0.0113 + 0.0006 


0.0151+0.0003 


0.0129 + 0.0001 



Table C.l: The transverse momentum distribution of the hard jet at the Tevatron. 



p T ,tt[GeV] 


H = m t /2 


dcr r fb l 
dp T ti L GeV J 

/i = m t 


IX = 2m t 


12.5 


16.43 + 0.07 


15.49 + 0.04 


13.2 + 0.2 


37.5 


37.25 + 0.05 


34.43 + 0.03 


28.46 + 0.02 


62.5 


13.56 + 0.02 


11.91+0.01 


9.61+0.01 


87.5 


5.66 + 0.01 


5.036 + 0.007 


4.053 + 0.004 


H2.5 


2.648+0.007 


2.373 + 0.005 


1.919 ±0.003 


137.5 


1.309 ±0.006 


1.1 94 ±0.003 


0.967 + 0.002 


162.5 


0.669 + 0.004 


0.623 + 0.002 


0.509 + 0.001 


187.5 


0.345 + 0.002 


0.332 + 0.001 


0.2720 + 0.0007 


212.5 


0.180 + 0.002 


0.1791+0.0008 


0.1482 + 0.0005 


237.5 


0.096 + 0.001 


0.0953 + 0.0006 


0.0797 + 0.0004 


262.5 


0.0466 + 0.0009 


0.0525 + 0.0005 


0.0443 + 0.0002 


287.5 


0.0243 + 0.0007 


0.0280 + 0.0003 


0.0239 + 0.0002 


312.5 


0.0113+0.0005 


0.0144 + 0.0003 


0.0126 + 0.0001 



Table C.2: The transverse momentum distribution of the total tt system at the Tevatron. 
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if 1 ,tL*- J ^' v J 


ij = m, /2 

LA. llt/^l — 


da r ft 1 

An 1 (" ,.\ J J 

a /?T t ^ e * 


H = 2wt 


12 5 


4 1 2 + 02 


3 58 +0 01 


2 910 + 007 


37.5 


in 96 + 03 


9 62 + 02 


7 78 + 02 


62 5 


1431 +0 03 


12 68 + 02 


10 35 +0 02 


87.5 


14 11 + 03 


12 63 + 02 


10 58 +0 02 


112.5 


11 ^7 + o 03 


10 66 + 02 


8 72 + 01 


1 37 5 


8 53 + 02 


7 94 + 01 


a 59 + O 01 


162 5 


5 75 +0 02 


5 48 +0 01 


4 570 + O 008 


187.5 


3 59 + 03 


3 59 + 01 


3 019 + 008 


212.5 


2 23 + 02 


2 24 + 01 


1 923+0 006 


237.5 


1.30 + 0.01 


1.391 ±0.008 


1.189 ±0.004 


262.5 


0.76 + 0.01 


0.837+0.006 


0.721+0.004 


287.5 


0.43 + 0.01 


0.488+0.004 


0.428 + 0.003 


312.5 


0.237 + 0.009 


0.286 + 0.004 


0.254 + 0.003 


337.5 


0.139 + 0.006 


0.167+0.003 


0.146 + 0.002 


362.5 


0.067 + 0.005 


0.088+0.003 


0.082 + 0.002 


387.5 


0.032 + 0.005 


0.051+0.002 


0.047 + 0.001 



Table C.3: The transverse momentum distribution of the top-quark at the Tevatron. 













\x = mJ2 


yu = m t 


\x = 2m t 


-3.8 


6.6 + 0.2 


6.3 + 0.1 


5.2 + 0.1 


-3.4 


16.0 + 0.3 


13.7+0.2 


11.2 + 0.1 


-3.0 


33.5 + 0.4 


30.1+0.3 


24.1+0.2 


-2.6 


70.1+0.8 


62.5+0.3 


50.9 + 0.2 


-2.2 


136.4 + 0.8 


123.3 + 0.5 


100.3 + 0.3 


-1.8 


237 + 1 


219.6 + 0.9 


188 + 9 


-1.4 


363 + 1 


338.1+0.8 


279.6 + 0.5 


-1.0 


483 + 1 


451.8 + 0.8 


376.2 + 0.6 


-0.6 


546 + 2 


519 + 1 


437 + 2 


-0.2 


573 + 2 


543.8 + 0.9 


454.5 + 0.6 


0.2 


573 + 2 


534.3 + 0.8 


446.2 + 0.6 


0.6 


542 + 1 


496.6 + 0.9 


409.1+0.6 


1.0 


467 + 1 


420.7 + 0.8 


341.4 + 0.6 


1.4 


354+1 


310+1 


248.6 + 0.5 


1.8 


229 + 1 


197.2 + 0.5 


165+9 


2.2 


130.5+0.7 


109.4 + 0.4 


86.4 + 0.3 


2.6 


65.9 + 0.5 


55.3 + 0.3 


43.0 + 0.2 


3.0 


32.4 + 0.4 


26.6 + 0.2 


20.9 + 0.2 


3.4 


15.2 + 0.3 


12.4 + 0.1 


9.7+0.1 


3.8 


6.7 + 0.2 


5.7 + 0.1 


4.42 + 0.08 



Table C.4: The pseudo-rapidity distribution of the top-quark at the Tevatron. 
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yt 


/u = m t /2 


yU = m t 


yU = 2m t 


-1.8 


9.8 + 0.3 


8.5 + 0.1 


6.9 + 0.1 


-1.4 


103.5 + 0.8 


96.3+0.4 


78.5 + 0.3 


-1.0 


388 + 1 


368.4 + 0.8 


314 + 9 


-0.6 


826 + 2 


777 + 2 


645.0 + 0.8 


-0.2 


1144 + 2 


1063 + 1 


888 + 2 


0.2 


1131+3 


1035 + 1 


856.6 + 0.8 


0.6 


799 + 2 


716 + 1 


581.9 + 0.8 


1.0 


379 + 1 


329.4 + 0.9 


270 + 9 


1.4 


100.9 + 0.6 


82.9 + 0.3 


63.5+0.2 


1.8 


10.1+0.1 


7.68 + 0.08 


5.46 + 0.07 



Table C.5: The rapidity distribution of the top-quark at the Tevatron. 







S [fb] 






/u = m t /2 


yU = m t 


yu = 2m t 


-3.8 


1.5+0.1 


1.26 + 0.09 


1.0 + 0.1 


-3.4 


10.8 + 0.4 


9.4 + 0.2 


7.5+0.2 


-3.0 


41.1+0.7 


35.8 + 0.4 


28.5 + 0.3 


-2.6 


100+1 


88.0 + 0.6 


70.3+0.4 


-2.2 


192 + 1 


168.2 + 0.7 


135.0 + 0.5 


-1.8 


295 + 1 


266 + 1 


215.5 + 0.6 


-1.4 


397 + 2 


357 + 1 


293.3 + 0.6 


-1.0 


456 + 2 


422 + 1 


348.4 + 0.6 


-0.6 


481+2 


447 + 1 


380 + 9 


-0.2 


471+2 


449 + 2 


374.4 + 0.7 


0.2 


475 + 2 


448 + 2 


374.9 + 0.7 


0.6 


477 + 2 


445 + 2 


380 + 9 


1.0 


457 + 2 


421 + 1 


348.8 + 0.5 


1.4 


396 + 2 


359 + 1 


292.8 + 0.5 


1.8 


294 + 1 


265 + 1 


214.9 + 0.5 


2.2 


191+2 


167.9 + 0.6 


135.3 + 0.5 


2.6 


101.1+0.9 


89.0 + 0.5 


70.1+0.4 


3.0 


41.8 + 0.6 


36.1+0.4 


28.7 + 0.3 


3.4 


11.3+0.3 


9.8 + 0.2 


7.7+0.2 


3.8 


1.5+0.1 


1.16 + 0.09 


0.92 + 0.07 



Table C.6: The rapidity distribution of the hard jet at the Tevatron. 
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dt ietrGeVl 

if 1 ,JclL VJV ' * J 


u — ni\ /2 


dcr r fb n 
A n T ■ ( L r lP v J 




25 


o 


o 


o 


75 


4201 + 11 


4045 + 5 


3553+4 


125 


1762 + 5 


1635 + 2 


1405 + 2 


175 


863 + 3 


802 + 2 


685 + 1 


225 


454 + 3 


428 + 1 


368 + 8 


275 


254 + 2 


242.7 + 0.7 


210.9 + 0.5 


325 


144 + 1 


142.1+0.5 


124.9 + 0.4 


375 


85.9 + 0.7 


86.8+0.4 


76.2 + 0.3 


425 


51.5 + 0.5 


53.8 + 0.3 


47.7+0.2 


475 


31.3 + 0.6 


34.1+0.2 


30.6 + 0.1 


525 


19.1+0.5 


21.8 + 0.3 


20.0 + 0.1 


575 


12.6 + 0.2 


14.7 + 0.2 


13.01+0.08 


625 


7.4 + 0.2 


9.4 + 0.1 


9.06 + 0.06 


675 


4.7+0.1 


6.40 + 0.09 


6.07 + 0.05 



Table C.7: The transverse momentum distribution of the hard jet at the LHC. 



Pr.tttGeV] 


yu = m t /2 


dcr r fb l 
d P7A L GevJ 

fi = m t 


fi = 2m t 


25 


2390 + 10 


1457 + 6 


920 + 2 


75 


1512 + 11 


2418 + 5 


2524 + 4 


125 


1915 + 5 


1722 + 3 


1460+1 


175 


938 + 3 


843 + 2 


714+1 


225 


496 + 2 


454+1 


383.9 + 0.7 


275 


274 + 1 


254.4 + 0.7 


217.9 + 0.5 


325 


157 + 1 


148.7 + 0.5 


128.2 + 0.4 


375 


88 + 1 


89.2 + 0.4 


77.7 + 0.3 


425 


54.5 + 0.8 


54.9 + 0.3 


48.3 + 0.2 


475 


31.6 + 0.6 


34.1+0.2 


30.7 + 0.2 


525 


18.8 + 0.6 


21.4 + 0.2 


19.8 + 0.1 


575 


11.5 + 0.2 


13.3 + 0.2 


13.07 + 0.09 


625 


6.6 + 0.2 


9.2 + 0.2 


8.50 + 0.08 


675 


3.7+0.1 


5.91+0.07 


5.72 + 0.06 



Table C.8: The transverse momentum distribution of the total tt system at the LHC. 
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PT.ttGeV] 


yu = m t /2 


dcr |- fb n 
dp Tt 1-GevJ 

yu = m t 


/i = 2m t 


25 


995 + 3 


881+2 


740 + 2 


75 


2152 + 5 


1930 + 4 


1639 + 2 


125 


2029 + 6 


1858 + 3 


1597 + 3 


175 


1338 + 11 


1276 + 2 


1115 + 2 


225 


723 + 7 


744 + 2 


662 + 2 


275 


368 + 2 


400 + 1 


365.8+0.9 


325 


178 + 2 


209.3 + 0.8 


196.4 + 0.5 


375 


83 + 1 


110.5+0.6 


107.5 + 0.6 


425 


37.9 + 0.7 


58.4 + 0.4 


60.1+0.5 


475 


14.9 + 0.6 


31.8+0.3 


33.8 + 0.2 


525 


4.5+0.4 


17.5+0.2 


19.6 + 0.1 


575 


0.1+0.3 


9.6 + 0.2 


11.4 + 0.2 


625 




5.4 + 0.1 


6.6 + 0.1 


675 




3.0 + 0.1 


4.21+0.07 



Table C.9: The transverse momentum distribution of the top-quark at the LHC. 











lit 


[X = m t /2 


yU = m t 


fj. = 2m t 


-5.7 


0.27+0.04 


0.21+0.01 


0.184 + 0.008 


-5.1 


0.85 + 0.04 


0.71+0.03 


0.58 + 0.01 


-4.5 


2.58+0.07 


2.27 + 0.03 


1.92 + 0.04 


-3.9 


7.5+0.1 


6.55 + 0.05 


5.35+0.05 


-3.3 


18.3+0.1 


16.3+0.1 


13.54 + 0.07 


-2.7 


36.4 + 0.2 


32.8+0.1 


27.86 + 0.07 


-2.1 


56.1+0.2 


52.2 + 0.1 


45.05 + 0.09 


-1.5 


69.6 + 0.2 


66.7+0.2 


58.3 + 0.1 


-0.9 


70.4 + 0.5 


69.7+0.1 


61.64 + 0.09 


-0.3 


66.6 + 0.2 


66.7+0.1 


59.28 + 0.09 


0.3 


67.0 + 0.2 


66.5+0.1 


59.39 + 0.09 


0.9 


71.4 + 0.3 


69.5+0.3 


61.57 + 0.09 


1.5 


69.6 + 0.2 


66.6 + 0.3 


58.1+0.1 


2.1 


56.3 + 0.3 


52.8+0.3 


45.10 + 0.08 


2.7 


36.0 + 0.2 


32.6 + 0.2 


27.80 + 0.08 


3.3 


18.3 + 0.1 


16.3 + 0.1 


13.43 + 0.08 


3.9 


7.8+0.2 


6.57 + 0.05 


5.43+0.04 


4.5 


2.6 + 0.1 


2.28 + 0.04 


1.82 + 0.03 


5.1 


0.83 + 0.04 


0.70 + 0.02 


0.61+0.01 


5.7 


0.25+0.02 


0.20 + 0.01 


0.171+0.009 



Table C.10: The pseudo-rapidity distribution of the top-quark at the LHC. 
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^[Pb] 

ay t Lr J 




Vt 


// = m t /9 


// = m t 
A* " t t 


= 2/Ht 


-3 9 


n 093 +0 008 


075 +0 006 

V./ . V./ / _!_ V./ . V./ \_7 \J 


q 053 +0 004 


-3.3 


2 39 + 06 


1.97 + 0.04 


1 55 +0 02 


-2.7 


12.9 + 0.1 


11 55 +0 09 


9 52 + 05 


-2.1 


35.7 + 0.2 


33.3 + 0.1 


28 35 +0 07 


-1 5 


66 7 + 2 


63 3 + 2 


55 1 +0 1 


-0 9 


96 5 + 3 

/U.J _!— V7 ■ ..7 


92 9 + 2 


81 5 + 1 


-0.3 


114.1+0.6 


111.1+0.2 


97.7 + 0.1 


0.3 


115.6 + 0.3 


110.9 + 0.2 


97.7 + 0.1 


0.9 


96.5 + 0.2 


92.9 + 0.2 


81.4 + 0.2 


1.5 


67.0 + 0.3 


63.5+0.2 


55.0 + 0.1 


2.1 


35.6 + 0.2 


33.1+0.1 


28.32 + 0.07 


2.7 


12.8 + 0.1 


11.45+0.07 


9.50 + 0.06 


3.3 


2.32 + 0.08 


1.95+0.03 


1.56 + 0.02 


3.9 


0.095+0.006 


0.071+0.005 


0.050 + 0.003 



Table C.l 1: The rapidity distribution of the top-quark at the LHC. 







S [pb] 




Jjet 


jj. = m t /2 


fi = m t 


yu = 2m t 


-5.1 


0.05+0.04 


0.18+0.02 


0.12 + 0.01 


-4.5 


2.2 + 0.1 


2.09 + 0.07 


1.88 + 0.04 


-3.9 


9.8 + 0.2 


9.8+0.1 


8.34 + 0.07 


-3.3 


23.0 + 0.2 


22.2 + 0.1 


19.38 + 0.08 


-2.7 


39.1+0.3 


37.1+0.1 


32.21+0.08 


-2.1 


53.3 + 0.5 


50.5 + 0.1 


44.01+0.09 


-1.5 


62.6 + 0.5 


59.9 + 0.2 


52.2 + 0.1 


-0.9 


68.3 + 0.6 


65.1+0.2 


56.8 + 0.1 


-0.3 


71.1+0.3 


67.2 + 0.1 


58.8 + 0.1 


0.3 


70.7 + 0.3 


67.0 + 0.3 


58.8 + 0.1 


0.9 


68.1+0.3 


65.4 + 0.3 


56.8 + 0.2 


1.5 


63.0 + 0.3 


60.1+0.1 


52.1+0.1 


2.1 


53.6 + 0.3 


50.5+0.1 


44.0 + 0.1 


2.7 


39.1+0.3 


36.9 + 0.1 


32.18 + 0.09 


3.3 


22.8 + 0.2 


22.3 + 0.1 


19.3 + 0.1 


3.9 


9.9 + 0.1 


10.2 + 0.4 


8.34 + 0.06 


4.5 


2.09 + 0.09 


2.2 + 0.1 


1.88 + 0.06 


5.1 


0.12 + 0.02 


0.13+0.02 


0.12 + 0.01 



Table C. 12: The rapidity distribution of the hard jet at the LHC. 
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